Valley splitting in Si quantum dots embedded in SiGe 
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We examine energy spectra of Si quantum dots embedded into Sio.75Geo.25 buffers using atomistic 
numerical calculations for dimensions relevant to qubit implementations. The valley degeneracy of 
the lowest orbital state is lifted and valley splitting fluctuates with monolayer frequency as a function 
of the dot thickness. For dot thicknesses < 6 nm valley splitting is found to be > 150 /xeV. Using 
the unique advantage of atomistic calculations we analyze the effect of buffer disorder on valley 
splitting. Disorder in the buffer leads to the suppression of valley splitting by a factor of 2.5, the 
splitting fluctuates with « 20/xeV for different disorder realizations. Through these simulations we 
can guide future experiments into regions of low device-to-device fluctuations. 



Understanding and design of silicon nanometer-scaled 
electronic devices has regained significant interest. This 
interest is sparked by the experimental progress that 
enabled the reproducible construction of geometries in 
which electrons are confined in three dimensions to 
length scales of a few nanometers and the potential ap- 
plications of this technology to ultra-scaled traditional 
CMOS devices. Emerging application of Si nanostruc- 
tures for qubit implementations due to long spin re- 
laxation times [l|, 0, 0] imposes additional stringent re- 
quirements on energy spectrum engineering, including 
the precise control of valley degeneracy. The six-fold 
valley degeneracy of bulk Si is reduced to two-fold de- 
generacy when electrons are confined to two dimensions 
(2D), such as at Si/Si02 interface in mainstream MOS- 
FETs. Already decades ago it was recognized that there 
is a small splitting between the two valleys in the low- 
est subband[4]. Recently, calculations predicted that val- 
ley splitting in narrow (few nm) SiGe/Si/SiGe quantum 
wells can be of the order of 10-100 meV and should fluc- 
tuate rapidly with the well thickness [H|, @, 0, H[ • However, 
experiments [§, LLO, LlJj produced valley splitting about 2 
orders of magnitude smaller than that prediction, which 
has been explained [l2| by the disorders of the Si/SiGe in- 
terface and in the SiGe buffer. The experiments [13j and 
theoretical methods indicated that additional spatial con- 
finement will minimize the role of interface disorder and 
increase valley splitting. In this paper we investigate the 
role of SiGe buffer disorder on valley splitting and answer 
the fundamental question of the size and controllability 
of valley splitting for relevant experimental structures. 

Three dimensional (3D) confinement of electrons can 
be achieved by various techniques. Electrostatic surface 
gating of 2D gas provides relatively weak and smooth 
spatial confinement potentials. In contrast, 3D confine- 
ment by Si/Si02 interface produces sharp potential with 
Coulomb energies approaching room temperature [3, EH 
and large valley splitting [la]. Recently, an alternative ap- 
proach to 3D confinement has been demonstrated with an 
advantage of lithographically denned epitaxial Si/SiGe 
interfaces using post-fabrication regrowth[17]. In this 



FIG. 1: Placement of Si (yellow) and Ge (blue) atoms in a) 
fully ordered (Si-Ge and Si-Si bonds) and b) partially ordered 
(Si-Ge, Ge-Ge and Si-Si bonds) 8-atom supercells. 



case spurious charging effects [HI, related to the traps in 
Si02 or unpassivated interface can be avoided, yet retain- 
ing sharp confining potential. We will simulate such de- 
fined Si nanostructures in SiGe buffers and explore sizes 
relevant for qubit implementations. Simulation capabil- 
ities to represent structures containing 10 million atoms 
explicitly enable the atomic representation of the dot, in- 
terfaces and the SiGe buffer. Atomistic simulations also 
present a unique opportunity to vary the amount of the 
buffer disorder in order to attain detailed understanding 
of the physics of valley splitting, including its magnitude 
and fluctuations. The valley splitting is primarily defined 
by the smallest dimension of the device and our conclu- 
sions are applicable to any Si nanostructure denned from 
SiGe/Si/SiGe quantum wells. 

Calculations of the energy spectrum are performed us- 
ing the NEMO-3D general purpose code, which repre- 
sents each atom in the domain explicitly. The theory 
undelying the tool and its relevant benchmarks are given 
m Refs [l9|,[20[. The structure is defined on the relaxed 
(001) Sio.75Geo.25 substrate and the Keating valence- 
force field model is used to adjust atomic positions to 
minimize the strain energy. Calculations of electronic 
structure are based on the 20 band sp 3 d 5 s* tight-binding 
model. The quantum dot was modeled as a lx x ly x lz 
rectangle grown on 37 nm-thick substrate and embedded 
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FIG. 2: Energy levels in l x x 20 nm xlO nm Si dot embedded 
into ordered Sio.75Geo.25 buffer. Energies are referenced to 
the valence band point. Inserts show spatial distribution 
of wavef unctions for the lowest levels. 
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FIG. 3: Valley splitting for the lowest orbital level as a func- 
tion of the dot size for ordered (black), partially ordered (red) 
and disordered (blue) Sio.75Geo.25 buffer. Bars indicate std 
deviation for each point. An example of valley splitting dis- 
tribution for 100 realizations of buffer disorder is shown in the 
histogram for l x — 25 nm, blue curve is the Gaussian fit. 



into 27 nm-thick Sio.75Geo.25 buffer, l z < l x ,l y , where z 
is along the growth direction. We investigated the influ- 
ence of the buffer thickness on electronic structure, there 
were no significant changes for substrates t s > 30 nm and 
buffers t& > 20 nm. 

For 25% Ge we can generate various placements of Ge 
atoms in the Sio.75Geo.25 buffer, with fully ordered con- 
taining only Si-Ge bonds, partially ordered containing 
single Ge-Ge bond per 8-atom supercell in a fixed po- 
sition, and disordered having random placement of Ge 
atoms retaining 25% composition, see schematic in Fig.[U 

We start with the analysis of energy levels and valley 
splitting in a dot embedded into a fully ordered buffer. 
Evolution of energy levels for a l x x 20 nmxlO nm dot 
is shown in Fig. [2] (the actual dot thickness l z = 9.85nm 
= 72 monolayers). All levels come in pairs, both levels 
in the pair having similar wavefunction envelopes (each 
level is also double spin-degenerate, which has been con- 
firmed by calculations and will be ignored for the rest of 
the paper). The 3D representations of the envelope wave- 
functions at 20% value are shown for the lowest 6 levels. 
The two lowest levels have similar s-type wavefunctions 
and represent the same orbital state with different val- 
ley number. The energy difference between them we call 
valley splitting AjJ. The next two levels have one node 
and belong to the next orbital state. For l x < 25 nm 
the p z -type state has lower energy than p x - and p y -type 
states due to the combination of sizes and effective mass 
anisotropy. The p x -type level has the highest sensitivity 
to l x , as expected, and for l x > 26 nm its energy becomes 
lower than that of the p z -type state. Energy separation 
between the ground and the first exited orbital states 
SE « 8 — 10 meV is large enough to restrict qubit Hilbert 



space to the lowest orbital state at low temperatures. 

Valley mixing results from superposition of two 
counter-propagating waves reflected from the opposite 
Si/SiGe heterointerfaces of the dot. The phase difference 
of the two waves depends on the details of the interface. 
The strength of the mixing depends on the amplitude of 
the wavefunctions at the interfaces, A v oc |x(^)| 2 , where 
x(h) is the value of the envelope of the electron wavefunc- 
tion at the dot boundary [7| . For p^-type and d z -type (top 
curve in Fig. [2j) states wavefunctions are pushed toward 
z-heterointerface and valley splitting for these state are 
significantly larger than for the ground and p x - or p^-type 
states. 

The most interesting question which can be uniquely 
studied by atomistic calculations is the role of buffer dis- 
order. In Fig. [3] valley splitting of the ground level is 
plotted for a l x x 20 nmxlO nm dot as a function of the 
dot size l x for fully ordered, partially ordered and com- 
pletely disordered buffer. For fully ordered buffer the 
valley splitting is ~ 0.5 meV, consistent with analyti- 
cal calculations. The value does not change significantly 
with the dot size, which confirms that valley splitting is 
primarily determined by the smallest dimension. For par- 
tially ordered buffers we see a reduction of A® by 10%, 
while for fully disordered buffer A^ is reduced 2.5 times 
to ~ 0.2 meV. To investigate fundamental reproducibil- 
ity of A® we performed calculations for 100 realizations 
of the buffer disorder for each point. The histogram of 
A^ for l x = 25 nm dot is plotted in the right frame. The 
distribution is Gaussian, with standard deviation of 9.4 
/ieV, which is ~ 5% of Aj. The bars on the main plot 
indicate standard deviation for other dot sizes. 

Inter- valley mixing is very sensitive to the smallest di- 
mension of the dot, l z , and fluctuates with a monolayer 
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FIG. 4: Valley splitting for the lowest orbital level of 25 nm 
x 20 nm x l z nm Si dot as a function of the dot thickness in 
monolayers (ML). l z is calculated using 1ML^ 0.13707 nm. 
Thin line connects points 1 ML apart, thick lines connect 
points 2 MLs apart (open dots for even MLs and solid dots 
for odd MLs). Bars indicate std deviations for different dis- 
order realizations. Stars show percent of the wavefunction 
penetrating into the buffer in z direction. In the inset A° for 
ordered and disordered buffers are plotted. Dashed line is AjJ 
obtained analytically for the 2D case. 



(ML) frequency A v oc cos(fco^)> where ko = 0.82(27r/a) 
is the center of the valleys and a is the lattice constant. 
Valley splitting as a function of l z with ML resolution 
is plotted in Fig. [4] (black line), and bars indicate stan- 
dard deviation for different disorder realizations. It has 
been noted that QWs with odd and even number of MLs 
belong to different symmetry classes 0- Indeed, if we 
connect A® for even and odd number of MLs we obtain 
two similar curves which fluctuate with a period of « 8 
MLs and are out-of-phase with each other. The value of 
Ay(l z ) f° r the dot embedded into a disordered buffer is 
reduced by a factor of 2.5, as shown in the inset. For 
comparison we also plot valley splitting calculated for 
the 2D QW using envelope function method [7[ (dashed 
line), which coincides with our calculations for the or- 
dered buffer. Saturation of valley splitting for large l z , 
compared with the analytical dependence, is due to 
an additional lateral confinement. In Fig.JHstars indicate 
percentage of the wavefunction |x(^)| 2 which penetrates 
the buffer above and below the dot, the envelope of A^ 
follows |x(^)| 2 as a function of l z . 

To summarize, we calculate energy levels and valley 
splitting for a small Si dot embedded in a disordered 
Sio.75Geo.25 buffer. We find that buffer disorder leads 
to the suppression of valley splitting by ~ 2.5 and actual 



values fluctuate with standard deviation of ~ 20 /xeV. At 
the same time disorder limits the lowest valley splitting, 
which can reach zero for a perfectly ordered buffer for 
some dot thicknesses, and dots with valley splitting > 
150 fieV can be predictably designed from narrow QW 
(lz < 6 nm). 
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